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ABSTRACT 



O 
H 

! Motivated by the desire to image exosolar planets, recent work by us and 

others has shown that high-contrast imaging can be achieved using specially 
shaped pupil masks. To date, the masks we have designed have been symmetric 
with respect to a cartesian coordinate system but were not rotationally invariant, 
thus requiring that one take multiple images at different angles of rotation about 
the central point in order to obtain high-contrast in all directions. In this paper, 
we present a new class of masks that have rotational symmetry and provide high- 
contrast in all directions with just one image. These masks provide the required 
1CT 10 level of contrast to within 4A/-D, and in some cases 3X/D, of the central 
point, which is deemed necessary for exo-solar planet finding/ imaging. They are 
also well-suited for use on ground-based telescopes, and perhaps NGST too, since 
they can accommodate central obstructions and associated support spiders. 



Subject headings: none supplied 
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1. Introduction 

With more than 100 extrasolar Jupiter-sized planets discovered in just the last decade, 
there is now great interest in discovering and characterizing Earthlike planets. To this end, 
NASA is planning to launch a space-based telescope, called the Terrestrial Planet Finder 
(TPF), sometime in the middle of the next decade. This telescope, which will ultimately 
be either an interferometer or a coronagraph, will be specifically designed for high-contrast 
imaging. Earlier studies (Brown et al. (2002)) indicate that a 4m class coronagraph ought 
to be able to discover about 50 extrasolar Earth-like planets if the telescope can provide 
contrast of 10~ 10 at a separation of 3X/D and that a 4 x 10m class telescope ought to be 
able to discover about 150 such planets if it can provide the same contrast at a separation 
of AX/D. 

One of the most promising design concepts for high-contrast imaging is to use pupil 
masks for diffraction control. So far, work in this direction (see Kasdin et al. (2002a); Brown 
et al. (2002); Spergel (2000); Kasdin et al. (2002b, 2003)) has focused on optimizing masks 
that are not rotationally symmetric and thus provide the desired contrast only in a narrow 
annular sector around a star, but at fairly high throughput. A full investigation of a given 
star then requires multiple images where the mask is rotated between the images so as to 
ultimately image all around the star. In this paper, we propose a new class of rotationally 
symmetric pupil masks, which are attractive because they do not require rotation to image 
around a star. Such masks consist of concentric rings. As we show, they can provide the 
desired contrast with a reasonable amount of light throughput. However, it is not clear a 
priori how to support a concentric-ring mask. While the simplest mechanism would be to 
lay the rings on a glass plate, even a tiny amount of scatter from the glass could destroy the 
high contrast. It is better to use an opaque structure for support, maintaining the purely 
binary nature of the mask. We show here that a certain type of spider support can be used 
without adversely affecting the contrast in the desired working range. 

The paper is organized as follows. In the next section, we show how to design high- 
contrast concentric-ring masks as a two-step optimization process. The first step is to solve a 
linear programming problem for an optimal apodization. The second step exploits the "bang- 
bang" nature of the optimal apodization to create a starting point for a local nonlinear search 
for an optimal mask. With masks in hand, Section 3 is devoted to considering the impact of 
support spiders on the point spread function and in particular on finding spiders that preserve 
the high-contrast region. Finally, we discuss how to design masks with central obstructions. 
Such masks, with their spiders, could be useful for telescopes with on-axis designs such as 
NGST and a variety of ground-based telescopes. For these, the mask's spiders could be used 
to block the spiders supporting the secondary. 
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2. Concentric- Ring Masks 



The image-plane electric field produced by an on-axis point source and an apodized 
aperture defined by a circularly-symmetric apodization function A(yJ x 2 + y 2 ) is given by 



and r(x,y) and 8(x,y) denote the polar coordinates associated with point (x,y). Here, and 
throughout the paper, x and y denote coordinates in the pupil plane measured in units of the 
aperature D and £ and ( denote angular (radian) deviation from on-axis measured in units 
of wavelength over aperture (X/ D) or, equivalently, physical distance in the image plane 
measured in units of focal- length times wavelength over aperture (fX/D). If the apodization 
function AQ takes only the values and 1, then the apodization can be realized as an 
aperture mask. 

For circularly-symmetric apodizations and masks, it is convenient to work in polar 
coordinates. To this end, let r and 9 denote polar coordinates in the pupil plane and let p 
and (f) denote the image plane coordinates: 





s 



where 



S = {(x,y): < r(x, y) < 1/2, 9(x, y) G [0, 2tt]} 



(2) 



x = r cos 9 
y = r sin 9 



c 



PCOS0 

psin (f>. 



(3) 



Hence, 



rp(cos 9 cos 4> + sin 9 sin 0) 
rp cos(# — (j)). 



(4) 
(5) 



The electric field in polar coordinates depends only on p and is given by 



1/2 2- 



E(p) 




(6) 
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where J denotes the 0-th order Bessel function of the first kind. Note that the mapping 
from apodization function A to electric field E is linear. Furthermore, the electric field in 
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the image plane is real- valued (because of symmetry) and its value at p = is the throughput 
of the apodization: 

1/2 

E(0) = 2tt / A(r)rdr. (8) 



The point spread function (psf) is the square of the electric field. The contrast require- 
ment is that the psf in the dark region be 10~ 10 of what it is at its center. Because the 
electric field is real-valued, it is convenient to express the contrast requirement in terms of 
it rather than the psf, resulting in a field requirement of ±10~ 5 . 

The apodized pupil that maximizes throughput subject to contrast constraints can be 
formulated as an infinite dimensional linear programming problem: 

maximize E(0) 

subject to -1(T 5 £(0) < E(p) < 1(T 5 £(0), p iwd < p < p owd , (9) 

0<A(r)<l, 0<r<l/2, 

where p; w d denotes a fixed inner working distance and p ow d a fixed outer working distance. 
Discretizing the sets of r's and p's and replacing the integrals with their Riemann sums, the 
problem is approximated by a finite dimensional linear programming problem that can be 
solved to a high level of precision (see, e.g., Vanderbei (2001)). 

The solution obtained for p iw d = 4 and p ow( j = 60 is shown in Figure 1. Note that the 
solution is of a bang-bang type. That is, the apodization function is mostly or 1 valued. 
This suggests looking for a mask that is about as good as this apodization. Such a mask can 
be found by solving the following nonlinear optimization problem. A mask consists of a set 
of concentric opaque rings, formulated in terms of the inner and outer radii of the openings 
between the rings: 

[ro,ri] first opening (10) 
[72,73] second opening (11) 
[74,75] third opening (12) 

; (13) 

b"2m-2,?"2 m -i] m-th opening (14) 



With this notation, the formula for E(p) given in (7) can be rewritten as a sum of integrals 
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over these openings: 

m-l ^ 

E(p) = 2tt^ / J (2nrp)rdr, (15) 

r2k 

m—1 

= ( r 2fc+i^i(2vrr 2fe+1 p) - r 2k J 1 (2nr 2 kp)) (16) 

Treating the ?Vs as variables and using this new expression for the electric field, the mask 
design problem becomes: 

m—1 

maximize tt ^ {r 2 2k+l - r\ k ) 

k=o (17) 
subject to -10" 5 £(0) < E(p) < 1()- 5 E(0), p iwd <p< p owd , 

< r < n < • • • < r 2m _! < 1/2. 

This problem is a nonconvex nonlinear optimization problem and hence the best hope for 
solving it in a reasonable amount of cpu time is to use a "local-search" method starting 
the search from a solution that is close to optimal. The bang-bang solution from the linear 
programming problem can be used to generate a starting solution. Indeed, the discrete 
solution to the linear programming problem can be used to find the inflection points of A 
which can be used as initial guesses for the r k s. The particular local-search method we used 
is the first-author's LOQO optimizer, which is described in Vanderbei (1999). Figures 2 and 
4 show optimal concentric-ring masks computed using an inner working distance of 4 and 
two different choices of outer working distances (40 and 60). 

There is great interest in decreasing the inner working distance as this rapidly increases 
the number of possible planets one could discover via TPF. Or, put another way, it allows one 
to study the same set of targets but with a smaller aperture instrument. Using p iw( j = 3.5, 
we first tried to solve the linear programming problem with p ow d = 10. There is no solution. 
We next tried p ow d = 7 and found a solution to both the linear and nonlinear optimization 
problems. The resulting mask is shown in Figure 6. 

Figure 7 shows a mask with p iw( j = 3 and p ow( j = 4.8. The linear program with p ow( j = 5.0 
is infeasible. We have even created a mask with p iw a = 2.5 and p ow d = 3.1. In the interest 
of space, we don't show this last mask. 

Of course, the obvious issue is ring support. How can the rings be supported without 
compromising the high-contrast area of the psf, without significantly reducing throughput, 
and without using glass which would introduce unacceptable scatter? The simplest answer is 
to use a "spider" similar to those used, for example, to support a secondary mirror. However, 
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spiders create diffraction spikes which can destroy the high-contrast. In the next section, we 
investigate the effect of spiders and show how to make them in such a way as to preserve 
high-contrast. 



3. Using Spiders for Support 

In this section, we study the impact that ring-connecting spiders have on the masks 
computed in the previous section. Normally, one thinks of a spider as a uniform-width 
support piece used to hang things. However, there is no specific need that these supports 
be of uniform width. In fact, since the masks will be made by electroforming technology, 
such a restriction to uniform width is completely unnecessary. Instead of constant-width 
spiders, we study constant-angular-width spiders. That is, each spider vane is a narrow 
sector, starting from a point at the center and flaring out with radius. Incorporating N such 
spiders uniformly spaced, the electric field, expressed in polar coordinates, is given by 

E{p,<j>) = JJ e' 2mrpcos{e -^ A{r)rdrd9 (18) 



where 



m—l 



A (r) = EWil( r ). ( 19 ) 



k=0 



S = {(x,y): < r(x, y) < 1/2, 6(x, y) E 0} , (20) 



r 2nn a 2n{n + 1) a 
2' N 2 



e=U ^^-^^--^ . (2D 

n=0 

a denotes the width of a vane in radians, the notation [a, b] denotes the interval on the real 
line from a to b, and the notation 1^] denotes the function that is one on the interval [a, b] 
and zero elsewhere. 

The integral in equation (18) can be expressed in terms of Bessel functions using the 
Jacobi- Anger expansion (see, e.g., Arfken and Weber (2000) p. 681): 



^IX COS t 



i m J m {xY mQ . (22) 
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Substituting into (18), we get: 

E(p,<f>) = I ' I ^i m J m {-2Tirp)e im{e - 4 ' ) A{r)rdrde (23) 

s m 

= [ $^"V m (-27rrp)e-^ ( f e ime de \ A{r)rdr (24) 
o m \e / 

The integral over 6 is easy to compute: 



27r(n + l) c 

JV-1 N 2 



y e ime ^ = J eim9dd ( 25 ) 

c "=0 2 7rn . Q 



JV T 2 

2vr - iVa m = 



-fsin(jiVf) m = jN,j^0 (26) 



otherwise. 



Substituting this result into (18), yields 

1/2 

£(p,0) = y J (-2nrp)(2n - Na)A(r)rdr (27) 



o 

1/2 



/ ^ jN JM-^rp)e- ilN<t, -sm(jN^)A(r)rdr. (28) 

Lastly, suppose that iV is even and use the fact that J_ m (x) = J m {— x) = (— l) m J m (x) to 
get the following expansion for the electric field: 

1/2 

E(p,4>) = {2-K-Na) J J (2nrp)A(r)rdr (29) 

o 

1/2 

oo ' 

/ JjN^rp) cos(jN((f) - n/2))- sm(jN^)A(r)rdr. (30) 
i-i ^ J 2 



The first term, involving the integral of Jo, is identical, up to a constant factor, to the 
formula for the electric field in the absence of spiders. Hence, the impact of the spiders is 
two-fold: (a) the original electric field is reduced by 1 — Na/ (2n) which is just the fraction 



- 8- 



of the open area that remains uncovered by the spiders, and (b) there are additional terms 
containing higher order Bessel functions. 

For large N, the effect of the higher-order Bessel terms becomes negligible for small p. 
Indeed, for x < ^jA(m + 1), 

(r/2\ m+1 

£ J »M 5 Uri^' (31) 

(which itself follows easily from the alternating Taylor series expansion of the m-th Bessel 
function: J m (x) = XXo(~ ^-Y(x/2) 2l+m /(l\(m+l)\)). From this we use the Schwarz inequality 
to estimate the magnitude of the effect of the spiders: 

1/2 



/ JjN^rp) cos{jN(<i)-7T/2))-sm(jN?-)A(r)rdr 

1 x 

j'=i 
1 00 1 

2^ 7W + I 



2^(jiV + l)!V 2 

for p owd < a/4(A^ + l)/ir. Since the last bound is dominated by e np o^ 2 /2, it follows from 
the dominated convergence theorem that this last bound tends to zero as N tends to infinity. 

The convergence to zero of the terms in the sum on j is very fast. In fact, if N is set 
large enough that max < I < xPowd Jn(x) < 10~ 5 , then the j = 1 term dominates the sum of 
all the higher-order terms and is itself dominated by the J term. Figure 9 shows a plot of 
^50, ^iooj and Ji 50 . These three Bessel functions first reach 10~ 5 at 35.2, 81.0, and 128.1, 
respectively. 

Figure 3 shows the psf for the mask designed with p iwd = 4 and p ow< j = 60 and three 
different spider configurations. The configuration consisting of 180 spiders each spanning 
0.0025 radians completely preserves the full high-contrast region. Note that this spider fills, 
in total, 180 x 0.0025 = 0.45 radians, which is 7.2% of the open area. 

Figure 5 shows the psf for the mask designed with p iwd = 4 and p owd = 40 and three 
different spider configurations. The configuration consisting of 120 spiders each spanning 
0.003 radians completely preserves the full high-contrast region. Note that this spider fills, 
in total, 120 x 0.003 = 0.36 radians, which is 5.7% of the open area. 

For the mask designed with p iwd = 3.5 and p owd = 7 (shown in Figure 6), 50 spider vanes 
suffice to preserve the high-contrast of the dark zone. The resulting psf is indistinguishable 
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from the one shown in Figure 6. If each spider spans 0.01 radians then they fill, in total, 
50 x 0.01 = 0.5 radians, which is 8.0% of the open area. 

For the mask designed with p iwd = 3.0 and p owd = 4.8 (shown in Figure 7), again 50 
spider vanes suffice to preserve the high-contrast of the dark zone and the resulting psf is 
indistinguishable from the one shown in the earlier figure. 

None of the masks presented have a central obstruction. Hence, there are many spiders 
all coming to a point at the center. This central area is probably the hardest part of the mask 
to manufacture. One probably needs to add a small bead at the center as an attachment 
point for the spiders. Using the same inner and outer working distances (3.0 and 4.8) as for 
the mask in Figure 7, we experimented with adding a constraint to the optimization problem 
that forces there to be a small central obstruction. In the nonlinear optimization model, the 
constraint simply amounts to fixing r to a specified small value. We tried several values and 
found that we could fix r to 0.01 (i.e., a 2% obstruction) and still find an optimal solution 
with essentially an unchanged throughput. 

It is a trivial matter to make designs with larger central obstructions to accomodate, 
say, a secondary mirror in an on-axis design. However, once the central obstruction becomes 
large, one loses the ability to get tight inner working distances. Figure 8 shows a design 
for a 31% central obstruction, an inner working distance of pi wd = 10 and an outer working 
distance of p owd = 40. With p iwd = 9, the linear optimization problem has a solution but not 
the nonlinear one. With p iwd = 8, not even the linear optimization problem has a solution. 

Some TPF concepts involve an elliptical pupil geometry since this might provide a means 
to achieve improved angular resolution in realizable rocket fairings. The designs presented 
in this paper are given in unitless variables. When re-unitizing, a different scale can be used 
for the x and y directions. In this way, these designs can be applied directly to elliptical 
pupils. Of course, the high-contrast region of the psf will also be elliptical with the short 
axis of the psf corresponding to the long axis of the pupil. 

4. Final Remarks 

There are two parts to TPF: discovery and characterization. Discovery refers to the 
simple act of looking for exosolar planets. Characterization refers to the process of learning 
as much as possible about specific planets after they have been discovered. The masks 
presented here are intended primarily for discovery since a single exposure with these masks 
can discover a planet in any orientation relative to its star. However, once a planet is found 
and its orientation is known, some of the asymmetric masks presented in previous papers will 
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be used for photometry and spectroscopy as they have higher single-exposure throughput. 

An alternative to pupil masks is to use a traditional coronagraph, which consists of an 
image plane mask followed by a Lyot stop in a reimaged pupil plane. Recently, Kuchner and 
Traub (2002) have developed band-limited image-plane masks that achieve the desired con- 
trast to within 3X/D. However, this approach suffers from sensitivity to pointing accuracy. 
Nonetheless, the approach is promising. In the future, we plan to consider combining pupil 
masks with image masks and Lyot stops to make a hybrid design that hopefully will provide 
a design achieving the desired contrast with an even smaller inner working distance. 

In this paper we have only considered scalar electric fields. We leave the important and 
more complex issue of how to treat vector electric fields, i.e. polarized light, to future work. 
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Fig. 1. — The optimal apodization for p iwd = 4 and p ow d = 60 and the associated psf. 
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Fig. 2. — Top. A concentric-ring mask designed to provide high-contrast, 10 10 , from X/D = 
4 to X/D = 60. Throughput is 0.134 = 17.00%. Bottom. The associated psf. 




Fig. 3. — Psf for mask in Figure 2 with spiders. Top Row. 100 spiders each spanning 0.005 
radians. Second Row. 150 spiders each spanning 0.003 radians. Third Row. 180 spiders each 
spanning 0.0025 radians. 
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Fig. 4. — Top. A concentric-ring mask designed to provide high-contrast, 10 10 , from X/D — 
4 to X/D = 40. Throughput is 0.131 = 16.62%. Bottom. The associated psf. 




Fig. 5. — Psf for mask in Figure 4 with spiders. Top Row. 50 spiders each spanning 0.01 
radians. Second Row. 100 spiders each spanning 0.005 radians. Third Row. 120 spiders each 
spanning 0.003 radians. 




Fig. 6. — Top. A concentric-ring mask designed to provide high-contrast, 10 10 , from X/D — 
3.5 to X/D = 7. Throughput is 0.184 = 23.5%. Bottom. The associated psf. 




Fig. 7. — Top. A concentric-ring mask designed to provide high-contrast, 10 10 , from X/D — 
3.0 to X/D = 4.8. Throughput is 0.251 = 31.9%. Bottom. The associated psf. 
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Fig. 8. — Top. A concentric-ring mask accomodating a large central obstruction. High 
contrast of 10~ 10 , from X/D = 10.0 to X/D = 40. Bottom. The associated psf. 




Fig. 9. — The Bessel functions J 50 , Jioo, and Ji 50 . They first reach 10 5 at 35.2, 81.0, and 
128.1, respectively. 



